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We discuss the production of gravity waves from the fragmentation of a supersymmetric conden- 
sate in the early universe. Supersymmetry predicts the existence of flat directions in the potential. 
At the end of inflation, the scalar fields develop large time-dependent vacuum expectation values 
along these flat directions. Under some general conditions, the scalar condensates undergo a frag- 
mentation into non-topological solitons, Q-balls. We study this process numerically and confirm the 
recent analytical calculations showing that it can produce gravity waves observable by Advanced 
Laser Interferometer Gravitational- Wave Observatory (LIGO), Laser Interferometer Space Antenna 
(LISA), and Big Bang Observer (BBO). The fragmentation can generate gravity waves with an 
amplitude as large as Q Gw h 2 ~ 10 1 and with a peak frequency ranging from mHz to 10 Hz, 
depending on the parameters. The discovery of such a relic gravitational background radiation can 
open a new window on the physics at the high scales, even if supersymmetry is broken well above 
the electroweak scale. 

o 

\ I. INTRODUCTION 

Supersymmetry is an appealing candidate for physics beyond the Standard Model. A generic feature of the scalar 
potential in supersymmetric generalizations of the Standard Model is the presence of flat directions parameterized by 
some gauge invariant combinations of squarks and sleptons. At the end of cosmological inflation, the formation of 
a scalar condensate along the flat directions can have a number of important consequences In particular, it can 
be responsible for the matter-antimatter asymmetry generated via Affleck-Dine (AD) mechanism dark matter in 
the form of Q-balls [!, 0, [H, H, 0] , squarks and sleptons driven inflation d, Q , and a curvaton [l(| • Note that for an 
inflation and a curvaton mechanism to succeed the flat direction must dominate the energy density of the universe 
. until they decay into Standard Model baryons. 
' In general, a supersymmetric condensate is unstable. An initially homogeneous condensate can break up into lumps 
Q^ , of the scalar field, called Q-balls [ll|, under some very generic conditions All phenomenologically acceptable 
supersymmetric generalizations of the Standard Model admit Q-balls [12], which can be stable, or can decay into 
fermions [11, OJ]. In many cases, the origin of the instability can be traced to running of the mass of the flat direction, 
due to logarithmic corrections |lj. For the squark directions, the leading order correction is usually negative due to 
the negative contribution of the gaugino loops When the condensate oscillates, the negative mass correction gives 
rise to an average negative pressure, which triggers the instability in the condensate [l5| . There are modes whi ch g row 
exponentially, and the condensate starts fragments into lumps, i.e., Q-balls. There are many analytical [1, EB, E3 EBl 
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and numerical [IjJ, l2fj, |2l|, |22j, |23| studies of Q-ball formation and their interactions [24 1 . 

It was recently pointed out that the process of fragmentation can serve as a source of gravity waves [25( with a 
detectable amplitude and with a peak frequency ranging from from 1 mHz to 10 Hz. This range of frequencies will be 
■ explored by a combination of upcoming detectors, such as the Laser Interferometer Gravitational Wave Observatory 
(LIGO) |2(1, t he Laser Interferometer Space Antenna (LISA) (27|, the Big bang Observer (BBO) [II] and the Einstein 
Telescope [29j. The spectrum of gravity waves is peaked near the longest wavelength, of the order of the fragmented 
region or of the Q-ball size (corrected for the red shift). We will discuss how the fragmentation of the condensate 
yields gravity waves and we will present both analytical and numerical results. 



II. FLAT DIRECTIONS AND Q-BALLS 



When the Standard Model is augmented by the scalar fields that carry baryon and lepton numbers, the non- 
topological solitons, or Q-balls, appear in the spectrum of such a theory At the end of inflation, large Q-balls, 
whose vacuum expectation values (VEV) are aligned with the flat direction, can form by fragmentation of the AD 
condensate 0, Q • The largest amplitude of gravity waves is attained when the flat direction condensate density is 
comparable to the total the energy density [25|. Usually this is not the case in AD baryogenesis, because the baryon 
and or lepton number carried by the condensate is constrained by the present-day baryon asymmetry of the universe. 
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However, some flat directions, namely those with B — L = 0, are not constrained because the net (B + L) asymmetry 
is destroyed by the electroweak sphalerons. The baryon and lepton number violating operators can contribute to the 
destruction of the Q-balls and can prevent them from dominating the energy density of the universe. 

The Q-ball with a global charge Q has the following properties. The scalar field inside the Q-ball has the form 

$>(x,t) = <f>(x) exp(iwi), (1) 

where <f> is real, and to ~ W3/2 ~ 0.1 — 10 TeV in gravity mediated supersymmetry breaking models. The global charge 
of the Q-ball is given by 

Q = uj J dxcj) 2 (x). (2) 

As long as the field has a time-dependent phase, it is associated with a non-zero global charge. 

In addition, there are flat directions with both B = and L = 0. This means that, in the {Re $, Im $} plane, the 
field undergoes radial motion without phase rotations. This can happen in the case of a flat direction inflaton [f| , or 
in the case of a flat direction curvaton [Io| . Due to the lack of a net charge, the fragmentation process generically 
leads to Q-balls and anti-Q-balls [U l30l. l3l|. which eventually decay. Part of the condensate energy goes directly 
into exciting the gauge bosons and gauginos, which eventually thermalize the universe. In what follows we will refer 
to rotations and radial oscillations of the flat directions in the {Re $, Im $} plane, depending on whether the flat 
direction carries a net global charge. 

The supersymmetric flat directions, by virtue of their couplings to the Standard Model fields, receive radiative 
corrections [l[. These corrections depend on the type of supersymmetry breaking. A typical potential contains soft 
supersymmetry breaking mass term and higher order non-renormalizable terms which arises by integrating out the 
heavy fields above the cut-off scale M 1 : 



V = m! /2 |4| : 



/ $<* \ |$|2d-2 



Here we have included the baryon and lepton number violating operators that are essential for AD baryogenesis and 
which play a role in decay of Q-balls, as discussed below. We consider M ~ Mpi ~ 2.4 x 10 18 GeV. It can be seen from 
the analyses of Ref. [32l. l33l |34| that most flat directions in MSSM are lifted by monomials of dimension 4. The soft 
supersymmetry breaking mass term, m, is proportional to the gravitino mass, i.e. m ~ m 3 / 2 ~ 0(100) GeV-O(l) TeV 
in gravity mediated supersymmetry breaking scenarios. The coefficient if is a parameter which depends on the flat 
direction, and the logarithmic contribution parameterizes the running of the flat direction potential. The value of K 
can be computed from the Renormalization Group (RG) equations, which, to one loop, give negative corrections due 
to gaugino loops [l[: 



K—f^> (4) 



where mi« is the gaugino mass and TO| is the slepton mass 2 . In the next section we will show that, whenever the mass 
of the homogeneous condensate receives a negative correction, the condensate undergoes fragmentation, an instability 
which leads to formation of Q-balls and, possibly, anti-Q-balls [l|, [3l| . 



III. AMPLIFICATION OF FLUCTUATIONS 



Let us consider the growth of the fluctuations in the flat direction condensate, which can be rotating ]J, [a| or 
oscillating [1,0, [Hj. The fluctuations in the field tend to grow when the average pressure is negative 0, [jj [l2 JTJEj 



1 For a gauge mediated case the potential along the flat direction is different, as discussed below. However, the energy density of the 
condensate is usually much lower than that in the gravity-mediated case, hence the fragmentation in the gauge-mediated case does not 
generate a comparable amount of gravity waves. 

2 A potential of the form of eq. {3}, with a negative value of K, can be obtained for a generic inflaton potential which has couplings 
to fermions and bosons, where the fermions belong to a larger representation than the bosons. The value of K is determined by the 
Yukawa interaction, h&ipip [3ll |; K ~ Ch? /Wit 2 , where C is the number of fermionic loops and h is the Yukawa coupling. 
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[ill |20| , which can be the case for K < in eq. ©. For field values </> < M, we find: 



F(0)^im| /2 ^(^y «^+«. (5) 

where we assume <C 1. The equation of state for a field rotating in such a potential is 

K \K\ 

P ^2^K P ^-— P > (6) 

where p and p is a pressure and energy density of the scalar field, respectively. Evidently, the negative value of 
K corresponds to the negative pressure, which signals the instability of the condensate. A linear perturbation 
analysis [5|, [31[ shows that the fluctuations grow exponentially if the following condition is satisfied (see Appendix 1): 



k 2 n 2 

a 



- ■+■ 2m 2 3/2 K ) < 0. (7) 



Clearly, the instability band exists for ne gati ve K, as expected from the negative pressure arguments The 
instability band, k, is in the range [!, [H, [l4ll9l l3~fl | 

°<^<%^^ 2m t/2\K\, (8) 

where a is the expansion factor of the universe. The most amplified mode lies in the middle of the band, and the 
maximum growth rate of the perturbations, see eq. (|33[) in Appendix 1, is determined by a ~ \K\m 3 j 2 /2 [3l| . 
The initial growth of perturbations can be described analytically in the linear regime by eqs. (|33I34[) with a(t) ~ 
\K\m 3 /2^t/2. When Scft/tfio ~ 0(1), the fluctuations become nonlinear. This is the time when the homogeneous 
condensate breaks down into Q-ball. We will study this process numerically. 

The energy density in the condensate depends on the model, and, foremost, on the type of supersymmetry breaking 
terms that lift the flat direction. This is because the potential along the flat direction depends on supersymmetry 
breaking (it vanishes in the limit of exact supersymmetry), and there are many ways to break supersymmetry. In the 
gauge-mediated supersymmetry breaking scenarios the potential can have the form [j| 

V(<l>)Km%]og(l + K) . (9) 



Here ms is the scale of supersymmetry breaking, which is of the order of 1 TeV. 

The main difference between gauge and gravity mediated cases for us is the mass per charge stored in the flat 
direction condensate and in the Q-balls that form eventually as a result of the fragmentation. In the gravity mediated 
scenarios, the mass density is po ~ m^ 2 cj>Q, the global charge density is uq ~ to 3 / 2 </>§, and the mass per unit global 
charge is of the order of m^, independent of the VEV <j) , for a review see In gauge mediated scenarios, where the 
Q-ball radius R ~ m^Q 1 ^ 4 , to <~ msQ 1 ^ 4 , 0o ~ msQ 1 ^ 4 , the mass density is p ~ m%, the global charge density is 
uq ~ w0q, and the mass per unit global charge is Po/uq ~ m|/0o [EEl- Furthermore, in a gauge mediated scenario 
the flat direction condensate never dominates the energy density of the universe. In what follows, we will concentrate 
on gravity mediated supersymmetry breaking scenarios. 



IV. GRAVITY WAVES 



The gravity waves are generated because the process of fragmentation involves inhomogeneous, non-spherical, 
anisotropic motions of the scalar condensate. As a result, the stress energy tensor receives anisotropic stress-energy 
contribution. The fragmentation of the condensate 3 takes place on spatial scales smaller than the Hubble radius. 
The gravity waves are generated at the time of fragmentation, i.e. when the linear perturbation in the flat direction 
condensate starts growing, as in eqs. (|33I34[) . 



3 Generation of gravity wave from the coherent oscillations of the inflaton field has been studied in Refs. \'Ml I'iTf . 
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In calculating energy density of the gravitational waves, we follow the transverse-traceless (TT) components of 
the stress-energy momentum tensor. By perturbing the Einstein's equation, we obtain the evolution of the tensor 
perturbations [3a |: 

V 2 

hij + 3Hh tJ - —h l3 = 167rGny , (10) 

where diHij = 11^ = and dihu = ha = 0. The TT part of the spatial components of a symmetric anisotropic 
stress-tensor can be found by using the spatial projection operators, Pij = 8ij — kikj, with fcj = fej/fc: 

Iljj (Aj) = A-ij,mn(^')^mn(A') 5 (H) 

where Ay,„ m (fc) = {Pi m (k)Pj n (k) — (l/2)i^j(fc)P mn (fc)V The TT perturbation is written as ft>y(t, k) = 
Ay,i m (fc)u i:) (i,k), where 

ilij + SHiiij — V 2 Uy = l&irGTij . (12) 

The source terms for the energy momentum tensor in our case are just the gradient terms of the flat direction 
condensate, 

T i:j = -^(V^iV^i + Vi&Vjfc) , (13) 

where 0i and 02 represent the real and imaginary parts of <fi, respectively. The gravitational wave (GW) energy 
density is given by 



PGW = 3^GV I d ^ h v h *iP ( 14 ) 



where V is the volume of the lattice. To estimate the magnitude of gravitational wave energy density on a lattice, we 
approximate eq. (|14[) by: 

PGW ^3^GV j^^uty (15) 
In numerical calculations we track the evolution of Haw using the u\2 and U21 components 4 . 

V. ANALYTICAL APPROXIMATIONS 

Let us estimate the fraction of energy density stored in the gravitational waves produced by the fragmentation of 
the condensate. The best-amplified mode is k 2 — fc^ aa . « m^l/CKl — |i^|/4) and the maximal growth rate of the 

fluctuation Scf> ~ Scf>oe a ^ +lkx is a — \K\m 3 / 2 /2. Therefore we obtain: 

{8<t>)=a5cj>=^\K\m m 5(t>, (16) 



During the fragmentation process the field perturbation grows from the initial value 

S(f>(t,x) = 54>o exp(m 3 / 2 \K\t/2 + ikx). 
Using the above approximation in eq. (|12p , we obtain 



2 

™3/2« - k liax U ~ k 2 max {5 (j){t , X )f . (17) 

1V1 Pl 



4 We have verified the approximate equivalence between Uu and hij by comparing the numerical results for the gravity waves obtained 
in the two approaches. 
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FIG. 1: Evolution of the condensate in two dimensions, calculated on the N = 64 lattice. White areas are the instability 

l 3/2- 



regions which eventually form Q-balls. The fragmentation occurs on the time scale t ~ O(10 2 - 10 a 




(a,)m 3/2 t = 900 



(b)m 3/2 t = 1050 



(c)m 3/2 t = 1200 



(d)m 3/2 t = 1350 



FIG. 2: Three-dimensional evolution of the condensate on the N = 64 lattice. Lumps of Q-matter are shown in color. The 
Q-matter breaks up into isolated lumps, Q-balls. 



Therefore, the fastest growing mode can be approximated by: 

u(t) a 2M p j 2 (6(t)(t,x)) 2 . 

The fragmentation becomes non-linear, and Q-balls form, when d<f>/(j> ~ 0{\). The majority of gravity waves are 
produced over the time interval determined by the condition; At ~ 2/(\K\) hx{(j> / 8(f)) . The final (saturated) gravity 
wave energy density can therefore be estimated by using eq. (|15p as: 



1 2 (duV \K\ 2 ml /2 4,(ty 



^~i M -u; — (18) 



The fractional energy density is then given by: 



Pgw , K , 2 ( 4>{t)\ 2 

\JW m \m7J • (19) 



i 2 

L 3/2 



Note that for a condensate, which is dominating the energy density of the universe at the time of fragmentation, the 
Hubble expansion rate is very small, i.e. m 3 / 2 ^> H(t), see [H, Q. For physically motivating parameters, we have 
chosen; m 3 / 2 ~ 100 GeV, <j>{t) ~ 10 16 GeV, and H(t) ~ 1 GeV, therefore, for a reasonable value of K ~ 0.1, we 
obtain, Ogw ~ 10 -6 . Note that View depends on the value of K, for K = there are no excitations of gravity waves. 
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VI. THE FATE OF Q-BALLS 

The formation of Q-balls is usually considered in the context of Affleck-Dine baryogenesis 0, 0, H[ • In this case, 
the energy density stored in the scalar condensate is small because it is related to the small baryon asymmetry of 
the universe. However, there is no reason why the supersymmetric scalar condensates could not have a much higher 
density if they carried a zero (B — L) charge, and if the Q-balls, formed via fragmentation, decayed before they came 
to dominate the energy density. The former requirement is sufficient for the primordial condensate to be independent 
from the baryon asymmetry of the universe because the net (B + L) global charge is erased by the sphalerons in the 
course of the electroweak phase transition. The latter has to do with the fact that, if the Q-balls forming from the 
fragmentation of the scalar condensate are long-lived, they can come to dominate the energy density in the universe 
causing an epoch of matter-dominated expansion that may not allow the efficient reheating at the end [38| ■ In the 
MSSM, there are flat directions that have B — L = 0, for example, QQQL, uude, QQud, QLue, etc. 

To estimate the range of the lifetimes of Q-balls, one must consider several decay modes. First, the scalar fields 
can evaporate into fermions carrying the same global quantum numbers. The decay of Q-balls via evaporation [39[, 
as well as Q-ball melting at finite temperature [7|, are both suppressed by the surface-to- volume ratio and can lead 
to some very long decay times [38l |40| . However, in the presence of baryon and lepton number violating operators, 
the decay may proceed much faster, because the Q-ball is only as stable as the U(l) symmetry is good. 

Let us first consider higher-dimensional operators suppressed by the scale M ~ Mpi. The supersymmetry preserving 
baryon number violating operators are given by F terms with dimensions larger than equal to 5 and D terms with 
dimension larger than 6. For example, the following baryon and lepton number violating operators can be written as 
F-terms: 

j^QiQjQkL t \e2 + —. 

where i,j,k,l are the generations (i ^ k). These interactions cause the B,L violation by A(B — L) = and 
A(B + L) = —2. A (B + L)-ball can lose its (B + L) charge and disintegrate via 2^3 processes that have cross 
section of the order of a ~ l/M 2 . The decay rate is given by: 

r~A^~™WlW#y-3 /2 , (2D 



C d —QiQjQkLt\ 8 2 + —UiEjUkdile* + h.c. , (20) 



Q dt v 11 \M_ 



i pi 



where ~ j 2 \K\ 2 m 3 /24>o ^ s the number density of Q-quanta inside the Q-ball. Although the flat direction undergoes 
fragmentation, but not all the number density of <j> field goes into forming a Q-ball. Here we have provided a 
conservative estimation, note that for K = 0, the Q-balls do not form. The factor 7 ~ 0.1 denotes the formation 
time scale of Q-balls, which is roughly given by t^ 1 ~ 71-^1^3/2, numerically one can see from Figs. (]1I2I) that 
m 3/2 t ~ 0(100) for \K\ ~ 0.1. The initial VEV of the flat direction is O ~ 10 16 GeV and M ~ M Pl . The Q-balls 
decay when T _1 is of the order of the Hubble time, we find that the Q-balls decay at temperature T ~ 10 5 GeV for 
mass m 3 / 2 ~ 10 2 GeV. 

However, the estimate of decay time in equation (I2ip is based on simplifying assumption of incoherent particle 
interactions inside the Q-ball, which may be inapplicable to scalars in a coherent state, such as Q-ball. It may be more 
appropriate to treat the Q-ball decay semi-classically, as discussed by Kawasaki et al. [4l| . Kawasaki et al. considered 

baryon number violating operators in eq. (|3|) that arise from supersymmetry breaking terms: Am 3 / 2 (yjjj + h.c.^ 

Such operators are essential for the Affleck-Dine baryogenesis to work. They found numerically that, for d = 4, the 
time scale of the Q-ball decay is (4l1 |: 



r ~ 10-° ( ^ ) m 3/2 , (22) 



This time scale is comparable to the Hubble time when the plasma temperature is T ~ 10 7 GeV for m 3 / 2 ~ 100 GeV 

and K ~ 0.1. The Q-ball formation occurs when the energy density is of the order of (l0 8 GeV) 4 , and the Hubble 
parameter is much larger than the decay width: H 3> T. Thus, the rate of baryon number violating processes is too 
slow to affect the Q-ball formation, but it is fast enough for the Q-balls to decay before they can dominate the energy 
density of the universe. 



VII. NUMERICAL RESULTS 



We solved the equations of motion for the field, eqs. ([55H4T]) in Appendix 2, numerically on a three dimensional 
cubic TV 3 lattice for N = 64, along with the evolution of Uij given in eq. (fT2"|) . For the purpose of illustration, 
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200 400 600 800 1000 1200 1400 

FIG. 3: Evolution of flaw, where the dark and light shaded dots correspond to the snapshots shown in Fig. [T]and Fig. [2] 
respectively. The comparison illustrates the agreement between two-dimensional and three-dimensional calculations. (Our final 
numerical results are based on the three-dimensional calculation.) 




FIG. 4: The panel on the left shows the spectrum of the growing perturbations in the flat direction condensate, and the panel 
on the right shows the gravity wave spectrum, as a function of time, on a N = 64 lattice. The overall normalization is arbitrary 
in both the cases. 



we have chosen the soft supersymmetry breaking mass TO3/2 = 10 2 GeV, which is motivated by gravity mediated 
supersymmetry breaking scenarios. We have considered a range of values for K: K E {0, —0.1, —0.2, —0.05, — 0.025} 5 . 
The initial VEV is taken (f>o = 10 16 GeV. The initial condition for the time-dependent phase was set randomly, as 
expected at the end of inflation. The initial small fluctuations of the condensate are set to be 54>/4> ~ 10~ 5 , as 
expected in inflationary cosmology. 

In Figs. [1] and [5] we show the fragmentation of the condensate. One can see the initial snapshot of a (nearly) 
homogeneous condensate in Fig. la. The process of fragmentation leads to islands with growing density contrast at 
later times, as shown in Fig. 1. Similar plots with different snapshots are shown in the 3D case, in Fig. 2, where one 
can see the fragmentation of the condensate into Q-balls. We can see how certain modes in the band of instability 
begin to grow leading to an eventual fragmentation of the condensate, and how these modes stop growing. At the 
same time, there is mode-mode mixing which excites higher-fc modes, and also affects the fragmentation process as 
shown in Figs. [1] and [2] The fragmentation does not happen isotropically, and coherent, macroscopic, non-spherical 
motions of the condensate create a quadrupole moment which leads to the creation of gravity waves. 

In Fig. [3] we show the evolution of the energy density stored in the gravity waves with respect to the critical energy 
density of the universe at the time of production. We plot flow as a function of m 3 / 2 t for 4> = 10 16 GeV along with 
2D slices and 3D isosurface plots for selected times, which correspond to the snapshots shown in Figs. [1] and [2] The 



5 The value of K = 



—0.2 is relatively large but we use this value for comparisons in studying the growth of Q.qw 
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(a) The final amplitude of the gravity waves does not depend on 
the initial perturbations. Here 4>0 = 10 16 GeV and 
m 3/2 = 100 GeV. 



(b)The final amplitude of gravity waves saturates for different 
values of K . However, for K = 0, there is no fragmentation and, 
therefore, there is no increase in the gravity wave amplitude. The 
value of | K determines the growth rate of the gravity waves. 



Here <j> = 10 16 GeV and m 



3/2 



100 GeV. 



FIG. 5: Effects of varying different parameters on the evolution of the gravitational wave fractional energy density Qgw- 



dark grey dots in Fig. [3] are results of the 2D calculation, and the light grey dots are results of the 3D calculation, at 
the times that correspond to the snapshots in Figs.[U andO 

At the onset of fragmentation all the scalar modes are excited; however, as discussed above, only a narrow band 
of the initial perturbation spectrum becomes unstable. The band of instability obtained from numerical calculations 
is shown in Fig. 0J It agrees with our analytical estimate, sec cq. ([5]). The gravitational waves spectrum is similar 
to that of the scalar perturbations, as shown in Figs. QJi, 2)d. However the gravity wave spectra are flatter. We have 
multiplied the gravity wave power spectrum by k 3 to highlight the growing mode. We have found that the peak of 
the gravity wave spectrum is determined by the longest wavelength which is of the order of the time scale when the 
Q-balls form, i.e. m 3 / 2 t ~ O(10 — 100), see Figs. [T]and[21 The growth eventually stops, once the fragmentation is 
near completion 6 . 

In Fig. O we show the temporal evolution of the gravitational wave density for different values of the initial field 
perturbations, different values of K and initial angular motion. In the perturbation plot, Fig. [5^, we find that the 
resulting gravitational wave energy density is independent of the size of the initial perturbation in the flat direction. 
The second plot, Fig. [5)d, shows different values of K determines how quickly the fragmentation proceeds, or the 
growth in the gravity wave density. We note that, when K is set to zero, there is no growing mode, the condensate 
never fragments, and no gravity waves are produced. 



VIII. OBSERVABLE SIGNAL 



After they are produced, the gravitational waves are decoupled from plasma. Let us estimate the fraction of the 
critical energy density p c stored in the gravity waves today: 



fW(to) = \f) I ) (23) 



1.67 x 10" 5 /100 
h 2 \~ 



5 / 1 nn\ V3 

n GW « io- 



With our current simulation we are unable to determine the exact peak frequency for the gravity waves. We would require a larger box 
size, i.e. N=256, to demonstrate this initial frequency, which is beyond the scope of the current paper. For the moment we will restrict 
ourselves with an approximate analytical discussion on an observable frequency, see section VIII. 
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where clq and Hq are the present values of the scale factor and the Hubble expansion rate, and a* and are the 
respective values at the time of the fragmentation. The estimate £Iqw ~ 10~ 6 obtained from our numerical results, is 
within one order of magnitude of the analytical estimates in eq. (fT9|) . LISA can detect the gravitational waves down 
to Vtcwh 2 ~ 1CP 11 at mHz frequencies. One can estimate the peak frequency of the gravitational radiation observed 
today, which is determined by the initial frequency, /* ~ u)k/2it. We obtain 

(24) 




where we have assumed that a* rs a r h (which also implies that one can neglect the expansion of the universe during the 
oscillations of the condensate). The numbers of relativistic degrees of freedom are g s ^ ~ 300, g s $ ~ 3.36. Here the 
subscript "rh" denotes the epoch of reheating and thermalization, while the subscript "0" refers to the present time. 
The typical frequency of the gravity waves is be determined by the size of the fragmented regions. This is roughly 
given by the scale at which fragmentation happens, i.e. m 3 f 2 tf ~ O(10 — 100). This can be seen in our numerical 
results shown in Fig. [T] and Fig. [2] For T r h ~ 1 TeV (the value of the reheat temperature is determined when the flat 
direction is responsible for reheating the universe [1, E3]), the frequency is of the order of 10 -3 — 10~ 2 Hz, which is 
in the right frequency range for LISA [27|. A higher temperature T r h ~ 100 TeV corresponds to the frequency range, 
10 — 100 Hz. Signals in both of these ranges will be accessible to BBO [28| and Einstein Telescope [2{j]. Since the 
supersymmetry breaking scale is related to the energy in the condensate, as well as the reheating temperature, future 
gravity wave experiments could be in a position to probe supersymmetry breaking scale above 100 TeV, which beyond 
the reach of the Large Hadron Collider (LHC). 

The spectrum of gravity waves produced by fragmentation is expected peak near the longest wavelength, of the 
order of the Q-ball size, which is ~ ((0.1 — 0.01)m 3 / 2 ) . The relatively narrow spectral width will help distinguish 
this signal from the gravity waves generated by inflation [43], which are expected to have an approximately scale- 
invariant spectrum (and a smaller amplitude). It may be difficult, however, to distinguish this signal from that 
generated by a first-order phase transition (44J, because both of these sources are expected to produce a relatively 
narrow spectrum determined by the Hubble parameter at the relevant time in the early universe. The gravity waves 
could also be generated in the electroweak-scale preheating. However, the amplitude of gravity waves expected in 



such a scenario |35| is considerably lower than that from the fragmentation of the supersymmetric condensate. LISA 
and BBO will be able to distinguish the gravity waves produced by fragmentation from those of point sources, such 
as merging black holes and neutron stars, which have specific "chirp" properties [45j . 



IX. CONCLUSIONS 



We have shown that the gravitational waves of observable amplitude could be produced when a homogeneous 
supersymmetric flat direction condensate fragmented into small lumps, i.e., Q-balls. The instability is a generic 
prediction for a supersymmetric flat direction in the early universe. The origin of the gravity waves is in the non- 
spherical anisotropic motions of the condensate that result from the growth of small initial perturbations. 

Gravity waves with the energy fraction as large as f^cw^o)^ 2 ~ 10~ n can be generated with a peak frequency 
ranging from mHz to 10 Hz, depending on the reheat temperature, which can vary in the range 1 — 100 TeV. The 
signal in the mHz frequency range can be detected by LISA, while a higher frequency 1 — 10 Hz is in the range 
of LIGO and BBO. The spectrum of gravity waves is different from many astrophysical sources identified by their 
chirp frequencies and from the gravity waves generated during inflation, which are stochastic waves with a scale 
invariant broad wavelength spectrum. However, a first-order phase transition in the early universe can produce a 
similar spectrum of gravitational radiation [44| • 

An observable amplitude of gravity waves can be produced by a class of fiat directions which carry no net (B — L) 
number and which are close to dominating the energy density of the universe at the time of fragmentation. Such flat 
directions are not responsible for generation of the baryon asymmetry via the Affleck-Dine scenario, although they 
may undergo a similar cosmological evolution as the flat directions discussed in connection with the matter-antimatter 
asymmetry. The (B + L) asymmetry generated by these flat directions is destroyed by the electroweak sphalerons. 
Although the identification of the origin of the signal may not be unambiguous, detection of these gravity waves by 
LISA, BBO, or EINSTEIN could open a window on supersymmetry in the early universe even if it is realized at a 
very high energy scale. 
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XI. APPENDIX 1 



Let us neglect any baryon number non-conservation and consider the homogeneous mode $ = (4>e ie ) / 'y/2 that obeys 
the classical equations of motion and and fluctuations about this classical solution, <j> — > + 6(f) and 9 — > 9 + 59. The 
equations of motion yield 0, [3l| 



■AH$- 9 2 4> + V'(4>) = 0, 
cf>9 + 3H(f>9-\-2<p9 = 0, 



(25) 
(26) 



for the homogeneous mode, and 



3H5<f> ~ 
ZH<f>60 + 2 



V 2 

- 9 2 5<p-^ r t 
a 

■ ■ ■ 6 
- 2- 



> + V"ty)84> = 0, 



V 2 

^59 = 0, 



(27) 
(28) 



for the fluctuations. Furthermore, 



V'M = ml /2 <f> 
V"((f>) 



1 + K + Kl0S {^p) 



m 



3/2 



1 + 3K + K log 



\2M 2 J 



(29) 
(30) 



Due to the conservation of the global U(l) charge in the physical volume, the solution has the property 9(f> 2 a? = const. 
If the energy density of the scalar field dominates the universe, the homogeneous part of the field evolves as 



9 2 {t) 



.-3/(2+*) 

<Pt 

ao J 

a 



l 3/2- 



(31) 
(32) 



To find the most amplified mode, we use the ansatz: 



'„ ( t \\ -3/(2+K) 
ao / 

m ) S9 e a ^+ ikx 
a 



SO 



(33) 
(34) 



If a is real and positive, these fluctuations grow exponentially, become nonlinear, and form Q balls [5j|. Substituting 
these into eqs. ([2~Tl) and (f2"8"|) , one obtains the following dispersion relation (3l| : 



g + 2m 2 3/2 Ka- 6K ^ 2 + K '> 



-2m 3/2 a-^ 2 + K Uo 



3K 
2 + K' 



-M 



2m 3 / 2 d 



a + a^ + — + 



k 2 

2 



3K 
2 + K 



(4 - 3K)H 2 



Ha 



= 0, 



(35) 
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where we have set ao = 1. If one neglects the cosmological expansion and assumes a <C a 2 for simplicity, one can 
reduce ea. (|35p to an approximate simplified equation: 



a 2 + ^ + 2m 2 3/2 K 



m 3 / 2 a 



. 9 k 



= 0. 



The perturbations grow exponentially if Red > 0, which yields 

k ( k <-> \ 
^U + 2ml /2 K)<0. 



(36) 



(37) 



XII. APPENDIX 2 



For our numerical analysis it is convenient to define dimensionless variables: 

ip = (p/m 3 / 2 , k = k/m 3 / 2 , r = m 3 / 2 t, £, = m 3/2 x , 

which we will use to study gravity waves during the non-linear fragmentation of the condensate. 
Writing ip = (tpi + itf2)/V2, we obtain the equations for the homogeneous mode: 



1 + K + K log 



2M 2 



<Pi = 0, 



(38) 



(39) 



where h = H/m 3 / 2 , i = 1, 2, and the prime denotes the derivative with respect to r. For the fluctuations, one obtains 

: 0, (40) 



d 2 d k 2 

-j-j + 3h— + ^ + Vij 
dr 2 dr a 2 



( Sipi 



where Vij denote the second derivatives with respect to <fii and ifj : 



Vu = l + K + K\og 



2M 2 



2K 



V w = V 2X = 2K 2 PlLP2 2 . 



(41) 



[1] K. Enqvist and A. Mazumdar, Phys. Rept. 380, 99 (2003); M. Dine and A. Kusenko, Rev. Mod. Phys. 76, 1 (2004). 
[2] I. Affleck and M. Dine, Nucl. Phys. B 249, 361 (1985). 

[3] A. Kusenko, V. Kuzmin, M. E. Shaposhnikov and P. G. Tinyakov, Phys. Rev. Lett. 80, 3185 (1998); J. Arafune, T. Yoshida, 

S. Nakamura and K. Ogure, Phys. Rev. D 62, 105013 (2000); Y. Takenaga et al. [Super-Kamiokande Collaboration], Phys. 

Lett. B 647, 18 (2007); S. Cecchini et al. [SLIM Collaboration], Eur. Phys. J. C 57, 525 (2008); I. M. Shoemaker and 

A. Kusenko, Phys. Rev. D 78, 075014 (2008). 
[4] A. Kusenko, L. Loveridge and M. Shaposhnikov, Phys. Rev. D 72, 025015 (2005); A. Kusenko, L. C. Loveridge and 

M. Shaposhnikov, JCAP 0508, 011 (2005). 

[5] A. Kusenko and M. E. Shaposhnikov, Phys. Lett. B 418, 46 (1998). 

[6] K. Enqvist and J. McDonald, Phys. Lett. B 440, 59 (1998) [arXiv:hep-ph /9807269 . 
[7] M. Laine and M. E. Shaposhnikov, Nucl. Phys. B 532, 376 (1998). 

[8] R. Allahverdi, K. Enqvist, J. Garcia-Bellido and A. Mazumdar, Phys. Rev. Lett. 97, 191 304 (200 6). R. A llahverdi, K. En- 
qvist, J. Garcia-Bellido, A. Jokinen and A. Mazumdar, JCAP 0706, 019 (2007) |arXiv:hep-ph/0610134| . 

[9] R. Allahverdi, A. Kusenko and A. Mazumdar, JCAP 0707, 018 (2007). 
[10] K. Enqvist, S. Kasuya and A. Mazumdar, Phys. Rev. Lett. 90, 091302 (2003). K. Enqvist, A. Jokinen, S. Kasuya and 
A. Mazumdar, Phys. Rev. D 68, 103507 (2003). R. Allahverdi, K. Enqvist, A. Jokinen and A. Mazumdar, JCAP 0610, 
007 (2006); 

[11] G. Rosen, J. Math. Phys. 9 (1968) 996; 9 (1968) 999; R. Friedberg, T. D. Lee and A. Sirlin, Phys. Rev. D 13, 2739 (1976); 
S. R. Coleman, Nucl. Phys. B 262, 263 (1985) [Erratum-ibid. B 269, 744 (1986)]; A. Kusenko, Phys. Lett. B 404, 285 
(1997). 



12 



[12] A. Kusenko, Phys. Lett. B 405, 108 (1997); Phys. Lett. B 406, 26 (1997). 

[13] G. R. Dvali, A. Kusenko and M. E. Shaposhnikov, Phys. Lett. B 417, 99 (1998). 

[14] K. Enqvist and J. McDonald, Phys. Lett. B 425, 309 (1998). K. Enqvist and J. McDonald, Nucl. Phys. B 538, 321 (1999). 
[15] K. Enqvist, A. Jokinen and J. McDonald, Phys. Lett. B 483, 191 (2000). 

[16] M. Axenides, S. Komineas, L. Perivolaropoulos and M. Floratos, Phys. Rev. D 61, 085006 (2000). 

[17] A. Pawl, Nucl. Phys. B 679, 231 (2004). 

[18] M. C. Johnson and M. Kamionkowski. larXiv:0805. 17481 [astro-ph], 

[19] S. Kasuya and M. Kawasaki, Phys. Rev. D 61, 041301 (2000); S. Kasuya and M. Kawasaki, Phys. Rev. D 62, 023512 

(2000); S. Kasuya and M. Kawasaki, Phys. Rev. Lett. 85, 2677 (2000). 
[20] T. Multamaki and I. Vilja, Nucl. Phys. B 574, 130 (2000). T. Multamaki and I. Vilja, Phys. Lett. B 535, 170 (2002). 

T. Multamaki, Phys. Lett. B 511, 92 (2001). 
[21] K. Enqvist, A. Jokinen, T. Multamaki and I. Vilja, Phys. Rev. D 63, 083501 (2001). 

[22] E. Palti, P. M. Sa ffin and E. J. Copeland, Phys. Rev. D 70, 083520 (2004). M. I. Tsumagari, E. J. Copeland and P. M. Saffin, 

larXiv:0805.3233l [hep-th] . 
[23] L. Campanelli and M. Ruggieri, Phys. Rev. D 77, 043504 (2008). 

[24] T. Multamaki and I. Vilja, Phys. Lett. B 482, 161 (2000). T. Multa maki and I. Vilja, Phys. Lett. B 484, 283 (2000). 
[25] A. Kusenko and A. Mazumda r, Phys. Rev. Lett. 101, 211301 (2008) larXiv:0807.4554l [astro-ph]]. 
[26] http://www.ligo.caltech.edu/ 
[27] http://lisa.nasa.gov/ 

[28] G. M. Harry, et al., Class. Quant. G ray. 23, 4887 (20 06) [Erratum -ibid. 23, 7361 (200 6)]; V. Corbin and N. J. Cornish, 

Class. Quant. Grav. 23, 2435 (2006); [http://universe.nasa.gov/new/program/bbo.htmll 
[29] |http: / /www.et-gw.eu / 1 

[30] K. Enqvist, S. Kasuya and A. Mazumdar, Phys. Rev. Lett. 89, 091301 (2002). 
[31] K. Enqvist, S. Kasuya and A. Mazumdar, Phys. Rev. D 66, 043505 (2002). 
[32] M. Dine, L. Randall and S. D. Thomas, Nucl. Phys. B 458, 291 (1996). 
[33] M. Dine, L. Randall and S. D. Thomas, Phys. Rev. Lett. 75, 398 (1995); 
[34] T. Gherghetta, C. F. Kolda and S. P. Martin, Nucl. Phys. B 468, 37 (1996). 

[35] J. Garcia-Bellido, D. G. Figueroa and A. Sastre, arXiv:0707.0839 [hep-ph]. J. Garcia-Bellido and D. G. Figueroa, 
larXiv:0801.4109 [gr-qc]. 

[36] J. F. Dufaux, A. Bergman, G. N. Felder, L. Kofman and J. P. Uzan, Phys. Rev. D 76, 123517 (2007). 

[37] R. Easther, J. T. Giblin, Jr., E. A. Lim, Phys.Rev.Lett. 99:221301,2007. R. Easther, E. A. Lim, JCAP bf 0604, 010, 

(2006). R. Easther, John T. Giblin, E. A. Lim, Phys.Rev. D77, 103519, (20 08). J. L. Cornou, E. Pajer and R. Sturani, 

Nucl. Phys. B 756, 16 (2006); A. Mazumdar and H. S toica, larXiv:0807.2570l [hep-th] . 
[38] S. Kasuya and F. Takahashi, JCAP 0711, 019 (2007) arXiv:0709.2634 [hep-ph]]. 
[39] A. G. Cohen, S. R. Coleman, H. Georgi and A. Manohar, Nucl. Phys. B 272, 301 (1986). 
[40] S. Kasuya and M. Kawasaki, Phys. Rev. D 64, 123515 (2001). 

[41] M. Kawasaki, K. Konya and F. Takahashi, Phys. Lett. B 619, 233 (2005) [arXiv:hep-ph/0504105| . 
[42] R. Allahverdi and A. Mazumdar, JCAP 0610, 008 (2006). 

[43] A. A. Starobinsky, JETP Lett. 30, 682 (1979) [Pisma Zh. Eksp. Teor. Fiz. 30, 719 (1979)]. S. Dodelson, W. H. Kinney and 
E. W. Kolb, Phys. Rev. D 56, 3207 (1997); M. Maggiore, Phys. Rept. 331, 283 (2000); T. L. Smith, M. Kamionkowski 
and A. Cooray, Phys. Rev. D 73, 023504 (2006): larXiv:0802.1530"1 [astro-ph], 

[44] M. Kamionkowski, A. Ko sowsky and M. S. Turn er, Phys. Rev. D 49, 2837 (1994); C. Grojean and G. Servant, Phys. 
Rev. D 75 (2007) 043507 [arXiv:hep-ph /0607107| . C. Caprini, R. Durrer and G. Se rvant, Phys. Rev . D 77 (2008) 124015 
arXiv:0711.2593 [astro-ph]]. C. Caprini, R. Durrer, T. Konstandin and G. Servant, arXiv:0901.1661 [astro-ph]. 

[45] B. J. Owen and B. S. Sathyaprakash, Phys. Rev. D 60, 022002 (1999). 




100 200 300 400 500 600 700 800 



m 3/2 t 



